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In Weyl materials the valence and conduction electron bands touch at an even number of isolated 
points in the Brillouin zone. In the vicinity of these points the electron dispersion is linear and may 
be described by the massless Dirac equation. This results in nontrivial topology of Berry connection 
curvature. One of its consequences is the existence of peculiar surface electron states whose Fermi 
surfaces form arcs connecting projections of the Weyl points onto the surface plane. Band bending 
near the boundary of the crystal also produces surface states. We show that in Weyl materials band 
bending near the crystal surface gives rise to spiral structure of energy surfaces of arc states. The 
corresponding Fermi surface has the shape of a spiral that winds about the projection of the Weyl 
point onto the surface plane. The direction of the winding is determined by the helicity of the Weyl 
point and the sign of the band bending potential. For close valleys arc state morphology may be 
understood in terms of avoided crossing of oppositely winding spirals. 

PACS numbers: 73.20.-r, 73.20.At, 71.10.Pm 


The possibility of existence of zero band gap semi¬ 
conductors with linear electron dispersion near the con¬ 
tact points of conduction and valence bands in the Bril¬ 
louin zone and the stability of this state with respect 
to Coulomb interactions has been investigated exten¬ 
sively beginning with the work of Abrikosov and Be- 
neslavskiipJ Recently there has been a renewal of theo- 
retical®3 an( j experimental research 1 ®!!] 0 n Weyl rne- 
terials, in which the electron dispersion near the point 
of crossing of the two non-spin-degenrate bands corre¬ 
sponds to that of relativistic massless (Weyl) fermions. 
The corresponding Berry connection curvature 2 ^ is diver¬ 
gence free with the exception of the band contact points, 
where monopole “charges” are located. By Nielsen and 
NinomiygP^ such Weyl point come in pairs of opposite 
helicity. In Weyl materials the Adler-Bell-Jackiw chiral 
anomaly,® 27 first discovered in particle physics, may be 
realized. This gives rise to giant anisotropic negative 

magnetoresistance.® 1 7 l 18 l 25 l 

Weyl materials possess surface stateipES whose Fermi 
surfaces have the shapes of arcs connecting the projec¬ 
tions of the Weyl points onto the crystal plane of the sur¬ 
face. Their origin is also rooted in the nontrivial topology 
of the Berry phase connection. 15 Arc states are close rel¬ 
atives of chiral edge states in graphene.®® The latter 
also owe their existence to the structure of Berry curva¬ 
ture in graphene,^, and their spectrum (as a function of 
momentum along the edge) consists of a line connecting 
the projections of the two Dirac points onto the edge. 
Evidence for arc states was recently observed by angle 
resolved photoemission spectroscopy (ARPES) in several 
Weyl materials. 19 23 

Another common mechanism of surface state genera¬ 
tion is band bending. It is caused by the difference be¬ 
tween electron environments at the surface and in the 
bulk. In previous theoretical considerations of arc states 
breaking of particle-hole symmetry due to band bend¬ 
ing near the surface was ignored (breaking of particle- 
hole symmetry at the graphene boundary leads to a fi¬ 
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FIG. 1. Sketch of the energy dispersion of the bulk states (in¬ 
terior of the Dirac cones) and the surface arc states (inclined 
plane segment) in the absence of band bending. 

nite dispersion of these, otherwise dispersionless, surface 
states.®®) The resulting energy dispersion of arc states 
is depicted in Fig. [Tj 

In this letter we study the effects of band-bending near 
the surface of Weyl materials. For conventional semi¬ 
conductors surface states caused by band bending are 
characterized by closed Fermi surfaces, see right panel in 
Fig- @ We show that in Weyl materials surface states 
caused by band bending hybridize with the arc states 
to form a single band with spiraling energy dispersion, 
shown in the left panel of Fig. [2] The corresponding 
Fermi surface is shaped as a spiral that winds about the 
projection of the Weyl points onto the boundary, see the 
left panel of Fig. |4j 

We consider a crystal boundary perpendicular to the 
2 -axis, as illustrated in Fig. [3j In the two band approxi¬ 
mation an electron state with a fixed crystalline momen¬ 
tum parallel to the boundary, p, may be described by a 
two-component pseudospinor, 
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FIG. 2. Left: the spiral energy surface of bound states 
for band bending potentials obtained from Eq. 0 with pa¬ 
rameters 7 = 4.5, zq = 0.01. For clarity only the part with 
d £ (—7 t/2 , 7 t/10 ) is shown. The surface terminates at the 
Dirac cones of particle-like and hole-like states of the contin¬ 
uum. Right: parabolic energy surfaces of bound states for a 
conventional semiconductor. 
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FIG. 3. Schematics of band bending potential U(z) and 
Dirac cones close to the boundary. 


Assuming for simplicity of presentation an isotropicpS 
dispersion near the Weyl points we write the effective 
low energy Hamiltonian in the form 

H = U(z) + vp ■ ct — ihvd z a z . ( 1 ) 

Here p = (p x ,p y ) and er are Pauli matrices, and the 
band bending is described by the potential U(z). The 
eigenstates with energy e obey the corresponding Dirac 
equation, 

e + ihvd z -U(z ) -vp \ f u p {z) \ __ 

-vp e - ihvd z — U(z) J v p (z) ) ’ 

( 2 ) 

where p = p x — ip y , and p = p x + ip y . 

In a general situation the projections of different Weyl 
points onto the plane of momentum parallel to the 
boundary do not coincide. In this case the boundary 
conditions do not couple different valleys and can be pa¬ 
rameterized by a single phase x, 


(v p - e lx u p ) 


= 0 . 

z=z 0 


( 3 ) 


The pseudospinor amplitudes obey second order differ¬ 
ential equations, i.e. 


<(*)■ 


e - U(z) 


hv 


+ i 


U'(z) 

hv 
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Wp(z) = 0. (4) 


Here the derivative with respect to z is denoted by prime. 
As follows from (|2|) the v amplitude may be obtained from 
the u amplitude by a differential operation 


pvp(z) = [e - U(z) + ihvd z ] u p (z) . (5) 


The momentum dispersion of states bound to the sur¬ 
face, e(p), is determined by discrete spectrum of the op¬ 
erator in Eq. ([!]) with the boundary condition Eq. (|3|. 
Although our conclusions regarding the spiral structure 
of the dispersion of surface states, e(p), are quite gen¬ 
eral, the specific shape of the dispersion depends on the 
details of the band-bending potential. 

The energy dispersion relation near the boundary of 
the continuum depends on the long (as compared to the 
lattice constant) distance behavior of the potential U(z). 
The corresponding macroscopic electric fields are created 
by the double layer charge distribution near the surface 
and depend on the screening properties of the Weyl ma¬ 
terial. In this respect one needs to distinguish between 
two cases: i) For zero band gap semiconductors at zero 
temperature, T = 0, linear screening is absent and U(z) 
is expected to decay as a power of z, ii) For zero band 
gap semiconductors at finite temperature or semimet- 
alsP^ U (z) decays exponentially at distances greater that 
the screening length (which is much longer than the lat¬ 
tice constant). In either of these cases the spatial depen¬ 
dence of U(z) at long distance may be analyzed using the 
Poisson’s equation, U"(z) = —47T e 2 n(z), where e is the 
electron charge and n(z) is the density of electrons. 

Semiconducting S3 case. In this case the thermo¬ 
dynamic electronic density of states (compressibility) 
dn/dp vanishes in the bulk for an undoped system. As a 
result, linear screening is absent, and the potential U(z ) 
exhibits a power law fall off with the distance 2 from the 
boundary. 

For small coupling constants, a = e 2 /hv, the band¬ 
bending region near the surface hosts many bound states. 
In this case the spatial dependence of the electron den¬ 
sity may be determined using the Thomas-Fermi ap¬ 
proximation. In the undoped crystal the electrochemi¬ 
cal potential is measured with respect to the Weyl point, 
i.e. U(z) + p{z) = 0. Due to the linear dispersion the 
electron density depends on the local potential U(z) as 
n(z) = —gU(z) 3 /3n 2 h 3 v 3 , where g is the number (even) 
of the Weyl points in the Brillouin zone. This yields the 
Thomas-Fermi equation in the form 


U"(z) 


4 ga 

37 x(hv) 2 


U 3 {z). 


( 6 ) 


The relevant solution is given by 


This corresponds to a vanish ing current normal to the 
boundary, j z = vip^a^ P 1 l 33 l 


U(z) = hv-, 
z 



( 7 ) 
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The short-distance behavior of the confining potential 
may not be understood within the low energy framework 
employed here and is accounted for by the boundary con¬ 
dition which should be imposed at distance of the 
order of the lattice constant. 

The second order differential equation Eq. Q with 
the potential U(z ) given by Eq. 0 can be reduced 
to the equation for the confluent hypergeometric func¬ 
tion. Introducing the dimensionless variable £ = 

2 \/|p | 2 — ( e / v ) 2 z/h we can write the general solution 
of Eq. Q in the form 

m (0 = C n exp(-C/ 2 )u;«), ( 8 ) 


were w(£) is expressed in terms of the confluent hyper¬ 
geometric function $(a, c as 

w(C) = $ (a, c; C) + Ci C 1_c $ (a - c + 1,2 - c; C) (9) 

with a = 7 (e/^v 2 \p \ 2 — e 2 — , and c = — 2 ry. 

For bound states the value of the constant C\ is de¬ 
termined from the condition that the wavefunction de¬ 
cay exponentially at C, —)> oo. Using the large distance 
asymptotic of $( 0 , c; (), we find 


a* T(o*) T(c) 
c(l — c) T(a) T(c*) ’ 


( 10 ) 


where T(x) is the gamma function. The spectrum of sur¬ 
face states is determined by substituting this form of u 
given by Eqs. ([8])-(10l into Eq. ^ and using the bound¬ 
ary condition, Eq. ([ 3 ]). Since the latter is imposed at 
short distances z = Zo h/ |p|, which corresponds to 
(0 1, the result of this substitution simplifies to 


i(<l>- 8 -%) r(fl*)r(c) 2z7ln( 2|p > j~ ci cos g) _ 

r (a)T(c*) 


( 11 ) 


Here a = y(tan 9 — i), and we parameterized the en¬ 
ergy of bound states, |e| < u|p|, by the angle 0 , so that 
e = n|p|sin0, with 9 G (— 7 r/ 2 , 7 t/ 2 ), and introduced the 
momentum azimuthal angle </>, p = |p|(cos<^, sin</>). 

Eq. (11) has a clear interpretation. Vanishing of the 
normal current j z = vip* a z i/j in a stationary state implies 
that the ratio v p (z)/u p (z) is a pure phase factor. The 
condition of exponential decay of the solution at z —> oo 
requires that v p (z)/u p (z) (z-^oo = e *w +e_ ^). The re¬ 
maining factor in the left hand side of Eq. © describes 
the dynamical phase accumulated during the motion be¬ 
tween the boundary and z —> oo. By the boundary condi¬ 
tion ^ the product of these phase factors must be equal 
to e lx in the right hand side. 

In the absence of the band bending potential the dy¬ 
namical phase vanishes. In this case bound states exist 
only in half the momentum plane, (j>~X £ (0, 7 r), in agree¬ 
ment with the treatments ignoring the band bending ef¬ 
fects near the crystal surface. At 7 7 ^ 0 the dynamical 
phase accumulated due to the presence of the band bend¬ 
ing potential allows the phase difference <fi — x to depart 



FIG. 4. Left: Fermi arc shap e for the undoped Weyl semicon¬ 
ductor determined from Eq. (111. The spiral makes infinite 
number of turns about the Weyl point. Right: Dashed lines 
- Fermi arc shapes for two nearby valleys of opposite lielic- 
ity for a square well potential, Eq. ( |12| ). The spirals make 
a finite number of turns (equal to the number of the bound 
states) about the Weyl points (red dots). Solid lines - avoided 
crossing of the two spirals due to finite valley mixing. 


from the interval (0,7r) and produces the spiraling struc¬ 
ture of the energy surface. By Levinson’s theorem®— this 
phase is given by 27 t times the number of bound states 
in the corresponding confining potential. As a result, the 
number of turns made by the spiraling energy surface 
around |p| = 0 is equal to the number of bound states 
caused by band bending. 


Due to the slow 1/z fall off, the potential in Eq. |7]) 
hosts an infinite number of shallow bound states. This 
is reflected in the logarithmic divergence of the phase on 
the left hand side of Eq. (11) at small parallel momenta, 
|p| —> 0. The shape of the Fermi surface of bound states 
in an undoped Weyl material is determined by setting 
9 = 0 (zero energy) in Eq. (11). The resulting spiral, 
which makes an infinite number of turns in the p x — p y 
plane, is plotted in the left panel of Fig. [4j 


The spectrum of bound states with nonzero energy 
is given by Eq. (11) with 9 / 0. For each momen¬ 
tum p = p| (cos <p, sin <p) it yields a discrete spectrum of 
bound states. This defines a surface in three dimensional 
space with coordinates p x , p v and e. In order to visualize 
this surface we can view Eq. © as the expression for 
the azimuthal angle ^asa function of energy ( 9 ) and the 
absolute value of momentum |p|. This gives a parametric 
representation of the dispersion relation of bound states. 
A three dimensional parametric plot of p x = |p|cos</>, 
p y = |p| sin 0 and e = v|p sinf? as a function of two pa¬ 
rameters, |p| > 0 and 9 £ that follows from 

Eq. (fTTj) is presented in the left panel of Fig. [2j 


Semimetallucase. For potentials U(z) that decay 
faster than 1/z into the bulk the number of bound states 
is finite. Therefore, the spiral Fermi surface makes a finite 
number of turns about the p = 0 point. As an illustration 
we consider the simplest example of such a potential, a 
square well of depth U 0 , U(z ) = —U 0 9(z 0 ~z). The treat¬ 
ment of this case as outlined above is straightforward and 
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leads to the following counterpart of Eq. (11), 


gl ( 0 +e _ f ) g+ (\p\cosd/h + iU) tan(qzo) = gix ^ 
q + (|p| cos 9/h — iU) tan(gzo) 

where U = Uq/Hv is measured with respect to wave vec¬ 


tor and q = yU 2 + 2C/|p| sin 9/h — (|p| cos 0/h) 2 , with 
the same parameterization as before. The shape of the 


Fermi surface can be obtained by setting 9 to 0 in Eq. (12) 
and is plotted in the right panel of Fig. [4] The number 
of windings it makes in the p x — p y plane is set by the 
number of bound states in the well, cxUzq. The charac¬ 
teristic momentum scale p s of the spiraling region is set 
by the depth of the well, p s « Uq/v. 

Note that the direction of the winding of the spiral¬ 
ing energy surface is opposite for the opposite helicity 
of the Weyl valley or opposite sign of the band bending 
potential. Thus for a known valley helicity the winding 
direction of the spiraling arc may be used to infer the 
sign of the band bending potential. 

Close valleys. The Fermi arcs join Weyl valleys with 
opposite helicity. If the momentum scale of the helical 
structure, p s ss Uq/v is smaller than the momentum dis¬ 
tance A p between the (projections on the terminating 
surface) Weyl points then the spiral structures near each 
valley bay be considered independently, as above. In the 
opposite case, p s < A p the interaction between the val¬ 
leys becomes important. In many of the Weyl materials 
Weyl valleys form nearby pairs with intervalley distance 
A p significantly smaller than the Brillouin zone size. In 
this case the effect of the intervalley interaction on the 
helical structure of arc states may be studied using the 
Dirac equation. Valley mixing of surface states arises 
from: i) valley coupling in the bulk, and ii) intervalley 
scattering from the crystal surface. The first effect is cap¬ 
tured by the four band Dirac Hamiltonian, so that Eq. 0 
may be replaced witlPI 7 jj = U(z)+VT 3 p-cr+ATia x with 
t being the Pauli matrices in the valley subspace. In this 
formulation the current operator becomes j = ev'^'T 3 tr'S. 
The intervalley scattering at the boundary is accounted 
for by generalizing the boundary condition 0 that cor¬ 
responds to the vanishing of the current j z normal to the 
boundary, 


Vp — Alup. 


(13) 


Here u p and v p are the two component (in the valley sub¬ 
space) generalizations of the spinor amplitudes u p and 
Up, and M is a 2 x 2 unitary matrix whose nondiagonal 
elements describe intervalley scattering from the bound¬ 


ary. In the absence thereof the matrix reduces to the 
diagonal form 0 , M ij =5 ij e i v. 

In the absence of valley mixing the arc states form two 
intersecting spirals winding in opposite directions (due to 
the opposite helicity of the Weyl valleys), shown in the 
solid lines in the right panel in Fig. [3J Arc state mor¬ 
phology at weak mixing may be understood in terms of 
avoided crossing of these spirals. It yields two alterna¬ 
tive morphologies of arc states Fermi surfaces. Type I 
all Fermi lines confined to the vicinity of the pair of the 
Weyl nodes (with one arc connecting them and, possi¬ 
bly one or more closed lines nearby). Type II - a pair 
of arcs emanating from the coupled Weyl nodes to other 
Weyl pairs. In the four band Dirac model both cases 
are realized depending on the choice of parameters in the 
Hamiltonian and the boundary condition (131. 


To illustrate the effects of valley mixing let us for sim¬ 
plicity neglect the valley coupling A in the bulk. In this 
case the ratio of the spinor amplitudes v^/u J p in valley j is 
a pure phase factor. For the potentials we considered, its 
value at the boundary, e la p, is given by the left hand side 
in Eqs. (11) and (12 pH. Thus, according to the boundary 


condition (|13|) the spectrum of arc states is determined 

= 0. Parameterizing 


5ije v 


- M, 


by the equation det 
the matrix M in Eq. (13) as M = e Ix e* T2e , where the 


angle 0 describes the valley mixing, and \ is the phase 
characterizing the intra-valley reflection, c.f. Eq. 0 , and 
using the phase factors e iQ p for the square well potential, 
Eq. (12) we get the energy dispersion of mixed valleys. 


At zero valley mixing, 0 = 0 the Fermi arcs are shown 
by the dashed lines in the right panel of Fig. [4] Avoided 
crossings due to small valley mixing results in type II 
morphology (solid lines), which is consistent with “tad¬ 
pole” arc shape observed in Refs. [HH23J 

In summary, we studied the interplay between the con¬ 
ventional, band bending mechanism of surface state for¬ 
mation, and the topological mechanism leading to for¬ 
mation of arc states. It results in a spiral structure of 
energy dispersion of surface states near the Weyl nodes. 
For close valleys, the arc state morphology may be un¬ 
derstood in terms of avoided crossing of the spirals of in¬ 
dividual valleys. Our findings significantly affect the den¬ 
sity and the dynamics of electrons in the surface metallic 
channel. 
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